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In the present paper we extend the multiparameter couphng constant metamorphosis, also known as the gener- 
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preserves integrability and a simple formula relating the solutions of the original system to those of the transformed 
- - , one are given. 
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The modern theory of integrable finite-dimensional dynamical systems concentrates mostly on Hamiltonian systems. 
For the latter, the presence of symplectic or Poisson structure combined with the Liouville theorem leads to a num- 
ber of important simplifications, in particular regarding the relationship among symmetries and integrals of motion. 
However, this also leads to certain restrictions because while considering various transformations of Hamiltonian 
C/^ ■ systems one naturally wants the transformed systems to be Hamiltonian too. Such transformations are of interest 
' for two reasons: they enable one to reduce the lists of integrable systems of certain form resulting from various 
'. classification procedures, and, sometimes, to relate new integrable systems to the known ones. 
I ^ i ' In the present paper we explore a somewhat surprising situation when the absence of need to preserve the Poisson 

■ or symplectic structure considerably simplifies things for general (non-Hamiltonian) dynamical systems including 
, ODEs. This is precisely the case for the multiparameter coupling constant metamorphosis, also known as the gener- 
[ alized Stackel transform, see [El [U \T7\ [3] and references therein. The said transform maps a set of integrals of motion 

■ in involution into another such set. In essence, it interchanges the parameters present in the integrals of motion with 
. the values of those integrals, thus producing a new set of integrals of motion. Note that the dynamical variables, i.e., 
] the phase space coordinates, are not affected by this transformation. 

■ The above recipe works provided the original set of integrals of motion depends on some parameters in a nontriv- 

■ ial fashion, and enables one to produce a number of interesting new examples or, conversely, to relate certain new 
i integrable systems to the known ones, cf. |17i [3]. However, in the Hamiltonian case there is a catch: the associated 

. ^ ', transformation for the equations of motion can only be written down upon restricting these equations to a common 
^ ■ level surface of integrals of motion generating the transform under study and turns out to be a rather nontrivial 



^ , reciprocal transformation, see [17] for details. We need this reciprocal transformation precisely because we want the 
- - ■ transformed equations of motion to originate from the transformed Hamiltonian through the Poisson structure which 
should be left intact. 

Below we extend the multiparameter generalized Stackel transform to general continuous dynamical systems, for 
which there is no Poisson structure to preserve. In this case it is possible to jettison the reciprocal transformation 
described above and define the associated transformation of equations of motion in a much simpler fashion. These 
results are summarized in Theorem [T] and Corollary [2] for general continuous dynamical systems and the overdeter- 
mined systems of first order PDEs, respectively, and in Corollary [3] for ODEs. What is more. Proposition [T] shows 
that the transformed system inherits existence of a Lax representation from the original one. The most important 
advantage of abandoning the reciprocal transformations is a very simple relationship among solutions of the original 
system and those of the transformed one, see Remark 1 below. 

As an aside, note that the requirement of presence of parameters in the dynamical systems under study and in their 
integrals of motion is not as restrictive as it may seem at first glance because the parameters can often be introduced by 
hand, e.g. through changes of dependent and independent variables, and often simple transformations like translation 
or rescaling of dependent variables enable one to produce interesting examples, see e.g. Examples 2, 3 and 5 below. 
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2 Coupling constant metamorphosis for general dynamical systems 



Consider an open domain M C IC" (K = M or C) and a dynamical system on M, 

dx^/fit = x"^, . . . , x", oi, . . . , flfe), a = l,...,n, (1) 

where are coordinates on M, and aj € IK are parameters. Thus, de facto we have a A:-parametric family of 
dynamical systems but for the ease of writing we shall refer to it below as if it were a single dynamical system. The 
same convention will apply to its symmetries, integrals of motion, Lax matrices, etc. We deliberately use the local 
setting instead of the more global one (vector fields on manifolds, etc.) because of the subsequent necessity to invoke 
the implicit function theorem which almost inevitably forces one to consider things locally, as explained below. In 
what follows all objects are tacitly assumed to be smooth in all of their arguments. 

To ([I]) we can naturally associate a vector field on M which depends on the parameters t, oi, . . . , Ofc, 

^ dx°' 

Recall that a (smooth) function / = f(t, x^, . . . , x", ai, . . . , a^) is an integral of motion (or a first integral) for the 
dynamical system ([TJ if we have 

df/dt + X{f) = 0. (2) 

Consider another dynamical system on M, 

dx^/fir = y"(t, x"^, . . . , x", ai, . . . , afc), a = l,...,n, (3) 
and the associated vector field on M, 

Recall that y is a symmetry for ([1]) if we have 

(Px°' (Px°' 



dtdr drdt ' ' ' ' 

where the derivatives are computed by virtue of ([I]) and ([3]), or equivalently, 

dY/dt + [X, Y] = 0. (4) 

Here and below [•, •] stands for the Lie bracket of vector fields (i.e., the usual commutator of differential operators) 
unless otherwise explicitly stated. 

Let ([1]) have k functionally independent integrals Ii, . . . , such that 

det{\m/daj\\ij=i^,„^k)^0- (5) 

Then by the implicit function theorem the equations 

Ii{x,t,ai, . . . ,ak) = bi, i = 1, . . . ,k, X e M, (6) 

where bi are constants, are (in general only locally, see the discussion in the beginning of the next section) uniquely 
solvable w.r.t. Oj. Denote the solution in question as follows: 

Ui = ii{x,t,bi,. . . ,bk), i = 1,. . . ,k, X € M. (7) 

The reason for this notation will become clear in a moment. 

If is a geometrical object on M (a function, a vector field, a tensor field, a differential form, etc.) which may 
depend on the parameters Oj, then K will stand for the geometrical object obtained from K by substituting /j for 
for all i = 1, . . . ,k. We shall write this as 

K = K\ J J . 
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Thus, for instance, X is a vector field obtained from X by substituting Jj for Oj for alH = 1, . . . , A:: 

and the associated dynamical system reads 

= X"(t,x\...,x",6i,...,6fc), a = l,...,n. (8) 

The only exception from the above notational convention is made for the functions !{ which are not related to 
in the fashion described above. The reason for this apparent discrepancy is that Ij turn out to be first integrals for 
X, see Theorem [1] below. 

Note that X and li depend on the parameters bi, . . . , ht, and that li obviously are functionally independent. 

In analogy with [T7], we shall refer to the procedure of passing from ([1]), X and Jj to ([8]), X and Ij as to the 
k-parameter generalized Stdckel transform (or as to the k-parameter coupling constant metamorphosis) generated by 
/i, . . . , Ifc. We shall also say that ([U is Stdckel- equivalent to ([5]). 

Just as in the Hamiltonian setting [17j, we have the following duality: when applied to ([8|), the A;-parameter 
generalized Stackel transform generated by Ii , . . . , 1^ brings us back to the integrals Ii , . . . , and the system ([T]) we 
have started with. 

As a final remark, note that the above construction admits the following geometric interpretation. 
Consider the extended phase space N = M X P, where P C K^, an open subset of K.^, is the space where our 
parameters live: P 3 a = {ai, . . . , a^)^; here and below the superscript T indicates the transposed matrix. 
The dynamical system ([1]) can be naturally extended to N upon setting 

dai/dt = 0, i = l,...,k. (9) 

In other words, the extended dynamical system ([T])+([9]) on is determined by the same vector field X as the original 
system ([1]) but X is now treated as a vector field on N. Conversely, the original system ([1]) is recovered from the 
extended one upon fixing the values of parameters a,, i = 1, . . . , fe. 

The extended system under study admits, in addition to k integrals of motion Ij, j = 1, . . . ,k, which are there by 
assumption, the 'obvious' integrals of motion Oj, j = 1, . . . ,k (for a moment, we ignore the possibility of existence of 
further integrals of motion for ([I])). Thus, foliates, at least locally, into common level surfaces of these 2k integrals. 
We shall denote this foliation by J^2k- 

Roughly speaking, the fc-parameter generalized Stackel transform, defined above, interchanges aj and bj as well 
as Ij and Ij, i.e., we choose a different way to parameterize the leaves of J^2k- 

The transform under study turns the extended system ([H)+([9]) into (f8|)+ (fT0|l . where (fTO|l reads 

dbi/dt = 0, i = l,...,k. (10) 

While ([I]) +([9]) and (f8l)+ (fT0l) coincide, at least locally (cf. the assumptions regarding the applicability of the implicit 
function theorem in the next section), on any given leaf of J-2ki they are different when considered on A^ as a whole. 

Our primary interest is, however, in the dynamical systems on the original phase space M rather than on the 
whole N, and this is where things become even more nontrivial: the transformed dynamical system ([8]) on M arises 
upon fixing the new parameters bj rather than the old ones aj, and ([8]) is a restriction onto M of the transformed 
extended dynamical system (j8])+ (ll0p on A^. 



3 Main results 

The above considerations suggest that the generalized Stackel transform should preserve a number of integrability 
attributes of the original system (e.g. symmetries, integrals of motion, etc.), and in the present section we state and 
prove the relevant results. 

Here and below we use the following blanket assumption. We suppose that the domain M and the ranges of values 
of time t and of the parameters aj and bj, j = 1, . . . ,k, are chosen so that the implicit function theorem ensures that 
the system ([6]) has a unique solution with respect to a,, i = 1, . . . , A;. In general this means making M and the ranges 
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in question sufficiently small because of the local nature of the implicit function theorem, but we deliberately do not 
fully spell out here the relevant conditions as there exists a number of situations when they would be too restrictive, 
e.g. when Ij are linear in all aj, j = 1, . . . ,k (or, more broadly, when we have explicit formulas for Ij and one can 
employ these to specify M and the relevant ranges, as it is the case for the majority of interesting examples). 

Theorem 1 Let the dynamical system flp have k functionally independent integrals Ii such that ^ holds. Define 
li and the transformed quantities like X as above. 
Then the following assertions hold: 

i) the functions li, i = 1, . . . ,k, are functionally independent integrals for and we have 

(11) 

, Jm o-re functionally 



dei[\dli/dbj\\i^j=i_k) /O; 

ii) if Ji, . . . , Jra is another set of integrals for such that all integrals Ii, . . . ,Ik, Ji, 
independent, then Ii, . . . , 1^, Ji, . . . , are functionally independent integrals for (0); 
Hi) if Yi, . . . ,Yr are linearly independent symmetries for ^\) such that 

Yp{Ij) = for all p = 1, . . . ,r and j = 1, . . . ,k, 

then Yi, . . . ,Yr are linearly independent symmetries for and 

Yp{Ij) = for all p = 1, . . . ,r and j = I, . . . ,k; 



(12) 



(13) 



iv) if under the assumptions of (ii) and (Hi) the symmetries Yi, . . . ,Ys, where s < r, span an involutive dis- 

s 

tribution, i.e., \Yp,Yq] = ^ Cpg(ai, . . . , Ofc, Ii, . . . , /fc, Ji, . . . , Jm)i^ for all p,q = l,...,s, then the symmetries 

9=1 

Yi, . . . ,Ys also span an involutive distribution, [Yp, Yg] = YI ($q^gj for all p,q = 1, . . . , s, where (?pq = c^q{Ii, . . . , Ik, 

9=1 

bi, . . . ,bk, Ji, . . . , Jm)- If Cp<j G lEC are constants (in particular, they do not depend on aj, j = 1, . . . , k), and thus Yg, 
g = 1, . . . , s, form a Lie algebra, then Yg, g = 1, . . . , s, form an isomorphic Lie algebra. 



Proof. Consider the identities 



bi = Ii{x,Ii, . . . ,/fc). 



l,...,k, X € M, 



(14) 



that follow from ([7]). Obviously, dbi/dt + X{bi) = 0, as bi are constants, and hence acting by d/dt + X on the 
left-hand side and the right-hand side of (jl4p yields 



. I § + !(/.) 



dt 



+ X{h) 



k 

E 

i=i 



dt 



+ X{I, 



ai=Ii,...,ak=Ik i-1 
k 

ai=Ii,...,ak=Ik j=l 



+ E(§ + x(/,.) 



dt 



I, 



doi 



,k. 



a\—I\,...,ak—Ik 



ai=Ii,...,ak=Ik 

Here we used the chain rule and the fact that are integrals for ([T|), and therefore by ([2]) we have 



dli 



0. 



Using ([5]) we readily find that dij/dt + X{ij) = 0, j = 1,. . . ,k, so Ij are indeed integrals for ([8]). 
Eg. pip follows from the implicit function theorem. The functional independence of Ij is immediate. This 
completes the proof of part (i) . 
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Next, Ij and Jg are integrals of motion for ([T]) by assumption, so dij/dt + X{Ij) = and dJs/dt + X{Js) = 0, 
and using the chain rule shows that Jg are integrals of motion for ([8]) as we have 



= U, s = i, . . . , m. 

ai=-flvi«fc=-ffe 



The functional independence of Ii, . . . , I^, Ji, . . . , easily follows from that of Ii, . . . , 1^, Ji, . . . , Jm, and thus part 
(ii) is also proven. 

In a similar fashion, further taking into account (jl2p and bearing in mind that dYq/dt+ [X, Yq\ = by assumption 
as Yq are symmetries for ([T]), we obtain 



= 0, g=l,...,r, (15) 

ai=J'i,...,aj.=/ft 



as desired. The linear independence of Yq readily follows from that of Yq. Let us stress that if we drop the condition 
(jl2p . the right-hand side of (|15p is in general no longer obliged to vanish, and hence the quantities Yi will no longer 
be symmetries for X, cf. Example 1 below. 

Finally, to prove (iv) we note that, in complete analogy with ()15p . we have 



|vi,i-,l=(K,rj+|(r,(y|^-r,(yf|) 



ai=/i,...,afc=/fc 
r 



s=l s=l 

and the result follows. □ 

Informally, Theorem [T] states that if Ii and Js are integrals and Yj are symmetries for ([TJ , and ([5]) and ()12p 
hold, then Ij and are integrals and Yj are symmetries for ([8]), the Lie algebra of symmetries Yj is 'essentially 
isomorphic' to that of Yj, and (llip and (I13p hold. By a slight abuse of terminology, it can be said that the 
multiparameter coupling constant metamorphosis (or the generalized Stackel transform) preserves integrals of motion 
and symmetries that respect the generators Ij of the transform in question. Proceeding in the spirit of the proof 
of Theorem [1] it can be shown that the multiparameter generalized Stackel transform preserves invariant curves 
and surfaces, Darboux multipliers, Jacobi multipliers and other similar structures. Thus, under certain technical 
assumptions the transformed system ^ inherits the integrability properties of ([T|). 

It is now appropriate to recall the definition of extended integrability due to Bogoyavlenskij [4]: 

Definition 1 ([4j) A dynamical system ([ip is integrable in the broad sense if it has m functionally independent 
integrals of motion Ji, . . . , J^, where n > m > 0, and n — TTi linearly independent commuting symmetries 1^, . . . , Y^—^ 
such that Yi(Jj) = for all i = 1, . . . ,n — m and all j = 1, . . . , m. 

Theorem [1] implies that the generalized Stackel transform preserves extended integrability. Namely, the following 
assertion holds. 

Corollary 1 Let (CP he integrable in the broad sense, with the integrals Jj and symmetries Yp as in Definition d 
Consider a k-parameter generalized Stackel transform generated by the integrals Ij which are functions of Jg, s = 
1, . . . , m, i.e., Ij = Ij{Ji, . . . , Jm), 3 = ^, ■ ■ ■ , k, and assume that ^ holds. 
Then the transformed system is again integrable in the broad sense. 

Proof. To prove this corollary it suffices to notice that we can construct from Ji , . . . , a new set of functionally 
independent integrals, say, Ig, s = 1, . . . ,m, for ([T]), so that Ij for s < k are precisely the generators of the Stackel 
transform in question. Then by Theorem [1] the transformed quantities Ig, s = 1, . . . , m, and Yj, j = 1, . . . ,n — m, 
meet the requirements of Definition [T] for ([8]), if so do Ig and Yj for ([T]), and the result follows. □ 
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Remark 1. Unlike the case of Hamiltonian dynamical systems considered in |17j . where the reciprocal trans- 
formation was involved, we have a very simple recipe for relating the solutions of ([T|) to those of ([8|). Namely, if 
_ ;=;a(^^ ai, . . . , a^), a = 1, . . . , n, is a solution for ([T]) then 

= H"(t,ai,...,afc)|,^,-^^^j^^ ,^^^j^,a = l,...,n, (16) 

is an implicit solution for ([8]). 
In particular, if the formulas 

= ai, . . . ,ak,Ci, . . . , Cn), a = 1, . . . ,n, 
where Ci, . . . , C„ are arbitrary constants, define a general solution for ([T]) then the formulas 

= ai, . . . , afc, Ci, . . . , Cn)\t=t,a^=h,...,a,=h ' a = 1> • • • > (17) 

define an implicit general solution for ([8]). Using similar considerations one can also readily find out how an implicit 
or parametric (general or particular) solution of ([T|) transforms into an implicit or parametric (general or particular) 
solution of ([8j). 

As a final remark, note that the multiparameter generalized Stackel transform preserves (the existence of) the 
Lax representations. Namely, it is readily verified that the following assertion holds. 

Proposition 1 Let ([ip admit a Lax representation of the form dL/dt = [M, L], where L and M are N x N matrices 
that depend on t, x^, . . . , x", oi, . . . , and on a spectral parameter \, and [,] stands here for the commutator of 
matrices. 

Then the transformed dynamical system (0) possesses a Lax representation of the form dL/dt = [M,L], where 
L= L{t,x^,...,x'',ai,...,ak,\)\a,=l,,...,a^=I^^ M= M(t, x^, . . . , x", oi, . . . , a^, A)|^^^j^^ ^^^^^-^ . 

4 Examples 

To illustrate the above results, we start with the following easy example. 

Example 1. Consider the one-component nonstationary dynamical system (i.e., a first-order ODE), 

dx/dt = ax/{{x -tf + ax), (18) 

see equation 1.4.3-2.15 in [T5] . 

For now let us work over C. Then (jl8p admits an integral / = Inx -|- a/{t — x), and hence is integrable in the 
broad sense. 

Setting oi = a and bi = b we readily find that L = (b — lnx){t — x) is an integral of motion for the transformed 
equation 

dx/dt = (6 — In x) / (x — t + 6 — In x) , 

which therefore is also integrable in the broad sense. 

As we have already noticed in Introduction, while many dynamical systems of interest do not involve parameters, 
we can often introduce the parameters 'by hand', e.g., through translation or rescaling of dependent variables. 

Example 2. Consider the following dynamical system from Example 2.22 of [7] 

dx/dt = -2x'^ + 2z, dy/dt = -3xy, dz/dt = 4xz - 2x(2x^ - 9y^). 

which has an integral of motion of the form z — x'^ + . 

Upon rescaling the variable y, y (a/3)^/^y, we obtain the system 

dx/dt = -2x'^ + 2z, dy/dt = -3xy, dz / dt = Axz - 2x{2x'^ - 3ay'^) (19) 

with an integral of motion 

/ = z — x^ + ay^. 



6 



Setting ai = a and bi = b and applying the general theory presented above we find that 

~ b + x'^ — z 

2 

y 

is an integral of motion for the transformed system 

dx/dt = -Ix"^ + 2z, dy/dt = -3xy, dz/dt = -2xz + 2x^ + Qbx. (20) 

We now see that in the transformed system the right-hand side of the third equation is independent of y, so we have 
a decoupled subsystem for x and z, 

dx/dt = -2x^ + 2z, dz/dt = -2xz + 2x^ + 66x, 

i.e., (]20p is, in a sense, indeed a somewhat simpler object than the original system (I19p . Moreover, (j20p admits a 
symmetry y = yd/dy. However, as Y does not preserve /, Y{I) = —21 ^ 0, this symmetry has no counterpart for 
the original system (I19p . 

On the other hand, ()20p is easily seen to have another integral I2 = x^/2 + z + 261n(y), and hence ()20p is integrable 
by quadratures (and integrable in the broad sense). 

Indeed, upon restriction onto the common level surface I = Ci and I2 = C2 the system (f20]l boils down to a 
single ODE, 

dy/dt = Tyv'6(Ci-6 + C72y2_261ny), 

which is obviously integrable by quadratures, and hence so is (pOj) . 

Now, by Corollary [T] the above implies that ()19p is also integrable in broad sense (and integrable by quadratures, 
as we can readily obtain the general solution for (jl9p from that of (|20p using (jl6p and (jl7p ). Note that the counterpart 
of /2 for (dH) reads 

I2 = x^l2 + z + 2(ay2 - + z) ln(y) 
Example 3. For a somewhat more elaborated example, consider system 9.25 from [9] 

du/dt = —uv^ + u + dv/dt = u^v — u — v, dw/dt = — v?' 

which has two integrals of motion, + + Inw^ and w(uv — 1). 

Upon rescaling u {a\)^l'^u^ v — > (ai)^/^?;, w — )■ w/a2-, we obtain the system 

du/dt = —aiuv'^ + u + V, dv/dt = oiu^w — u — w, dw/dt = aia2{v'^ — u^) 

with the integrals of motion 

Ii = ai(u^ + u;^) + In(w^), /2 = 021^(01^^; — 1). 
Consider the two-parametric generalized Stackel transform generated by Ii and I2. We find that 

= (61 - \n{w'^))/{u^ + ly^), /2 = -b2{u^ + w"^) / {uvw{\n{w'^) - bi) + (n^ + i(j^)u;) 

are integrals of motion for the system 

du/dt = -{bi-\n.{w'^))uv'^/{u'^ + w'^) + u + v, 
dv/dt = {bi -\n.{w'^))u^v/{u^ + w"^) - u - V, 

dw/dt = b2{bi -\n.{w'^)){u^ - v'^)/{uvw{\n.{w'^) -bi) + {u^ + w'^)w). 
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5 Generalized Stackel transform for overdetermined partial differential sys- 
tems 

Theorem [1] admits a natural generalization to the overdetermined systems of first-order PDEs which naturally arise 
e.g. in the study of zero-curvature representations, pseudopotentials and Backlund transformations for integrable 
(systems of) PDEs, cf. e.g. [Ij and references therein. 

Namely, consider an overdetermined system of first-order PDEs of the form 

^ =X2(t\...,t'^,x\...,x",ai,...,afe), a = l,...,n, A=l,...,d, (21) 
and assume that this system is in involution, i.e., 

a = l,...,n, A,B = l,...,d, (22) 



where the derivatives are computed by virtue of (j2ip . or equivalently, 

dXA/dt^ -dXB/dt^ -[Xa,Xb] = 0, A,B = l,...,d. (23) 

Corollary 2 Let i f^) hold, and let (E2P have k joint (i.e., such that dlj/dt^ + X^(/j) = for all A and j) 
functionally independent integrals Ii, . . . ,Ik such that (0) is satisfied. 
Then the following assertions hold: 

i) the vector fields Xa, A = 1, . . . , again commute: [Xb,Xb\ =0, a,b = 1, . . . ,d, and hence the transformed 
system 

^ =X^(t\...,t'^,a;\...,x",6i,...,6fc), a = l,...,n, A = l,...,d, 
is again in involution; 

ii) the functions Ii, i = 1, . . . ,k, are functionally independent joint integrals for the vector fields Xa, a = 1, . . . ,d, 
and we have 

det(ldii/dbj\\ij=i_k) /O; 

Hi) if Ji, . . . , Jm is another set of joint integrals for Xa, A = 1, . . . ,d, such that all integrals /i, . . . , I^, Ji, . . . , 
Jm are functionally independent, then Ii, . . . , 1^, Ji, . . . , Jm are joint functionally independent integrals for Xa, A = 
l,...,d; 

iv) ifYi,...,Yr are linearly independent joint (i.e., dYq/dt^ + [Xa, Yq] = for all A and q) symmetries for Xa, 
A = 1, . . . ,d, such that 

Yp{Ij) = for all p = 1, . . . , r and j = 1, . . . , k, 
then Yi, . . . ,Yr are linearly independent joint symmetries for Xa, A = 1, . . . ,d, and 

Yp{Ij) = for all p = 1, . . . , r and j = 1, . . . , k; 

v) if under the assumptions of (Hi) and (iv) the symmetries Yi, . . . , Yg, where s < r, span an involutive distribu- 

s 

tion, i.e., \Yp, Yq\ = ^ Cpq(ai, . . . , ak,Ii, . . . , /fc, Ji, . . . , Jm)Yg for all p,q = 1, . . . , s, then the symmetries Yi, . . . ,Ys 
also span an involutive distribution, i.e., [1^,1"^] = Cpgl^, for all p,q = l,...,s, where (%q = Cpg(/i, . . . , Ifc; 

9=1 

5i, . . . , Ji, . . . , Jm)- If Cpq £ IK are constants (in particular, they do not depend on aj, j = 1, . . . ,k), and thus Yg, 
g = 1, . . . , s, form a Lie algebra, then Yg, g = 1, . . . ,s, form an isomorphic Lie algebra. 

Stating the counterpart of Proposition [1] for ([2T]) is left as an exercise for the reader. 
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6 Applications to ODEs 

Consider an ODE resolved with respect to the highest-order derivative: 

d"'u/dz"' = F {z, u, du/dz, . . . , d^'-^/dz"'-^ ,ai,...,ak). 

Let n = m, and put 



(24) 
(25) 



x^ = u, x^ = du/dz, x"" = d"'~\/dz"'-\ 

Consider a dynamical system 

dx^/dt = x^, dx^/dt = x^, ...,dx"'^^/dt = x"', dx"^ /dt = f{t,x\ . . . ,x"',ai, . . . ,ak); (26) 



here f{t,x^, . . . ,x™,ai, . . . ,ak) = F i z,u, 



du 



dz''"' dz'^-^ 



, ai, . . . , Ofc 



z=t,u=x^ ,du/dz=x^ ,...,d™- ^u/dz 



It is well known that the dynamical system (j26|) is equivalent to (|24p . and we can readily apply the result of Theorem [T] 
to ([26]). 

What is more, it is immediate that upon applying the multiparameter generalized Stiickel transform to ()26p we 
obtain the system of the same kind, that is, 



dx^/dt = x^, dx'^/dt = x^, . . . , 
dx'^-^/dt = X™, dx'^/dt = fit, x\ . . . , X™, 6i 
which is, through (j25p . equivalent to an ODE of the form 

d™u ^ / du d'^-^u 



1 • • • 1 ^K, J "i 



h), 



(27) 



dz'' 



du 
dz ' 



F { z,u,—,..., j_^_-^^ ,bi,...,bk 



dz' 



where 



' ■ ■ ■ ' Ij^m-^ '^i'-- -^^k) = fit,x\.. . . .,bk) 



(28) 



(29) 



Thus, we have obtained a transformation relating the ODEs (j24p and (j28p . and this transformation preserves the 
integr ability properties. 

In view of the particular interest in the study of ODEs let us restate Theorem [1] for this special case directly in 
terms of ODEs. To this end we first recall the relevant definitions following |14j . 

A generalized vector field Y = h{z,u, du/dz, . . . ,d^~^u/ dz'^~^)d / du is a (generalized) symmetry for (j24p if we 
have 



m— 1 



dF 



j=0 



Here uq = u,Uj = d^u/dz^ , and we have introduced the so-called operator of the total z-derivative 

D- — F ^ V — 

dz dum~i ^ ^^^duj 

(here we treat z and Uj as formally independent entities, see e.g. jTlj for details). 

Also, a function / = I{z,u,du/dz, . . . ,d^~^u/dz^~^) is a (first) integral for if D{f) = 0. 

It is easily seen that upon passing from (j24p from (j26p an integral of motion I and the prolongation (see e.g. 
of a symmetry Y 

m— 1 



pr 



Y=Y,D^{h) 



j=0 



become respectively an integral and a symmetry for (j24p in the sense of the definitions from Section [2j 
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If Yi = hi{z,u, du/dz, . . . ,d"^ ^u/dz"^ ^)d/du, i = 1,2, are two symmetries for ()24p in the sense of the above 
definition, their commutator is [14j given by the formula 

[Yi,Y2] = {piYi{h2)-wY2{hi))d/du,. (30) 

and of course it is again a symmetry for ()24p . 

With ah this in mind we are ready to state the ODE version of Theorem [TJ 

Corollary 3 Under the above assumptions, let \24^ he an ODE whose right-hand side depends on k parameters 
ai, . . . ,ak, and let [24^ have k functionally independent integrals h, . . . ,Ik such that ^ holds. 
Then the following claims hold: 

i) the functions li, i = 1, . . . ,k, are functionally independent integrals for the transformed ODE (28\) . and we have 

det(ldii/dbj\\ij=i_k) /O; 

a) if Ji, . . . , Jm is another set of integrals for [24\ ) such that all integrals Ii, . . . ,Ik, Ji, . . . , Jm are functionally 
independent, then /i, . . . , /^., Ji, . . . , Jm are functionally independent integrals for I128\} : 
Hi) if Yi, . . . ,Yr are linearly independent generalized symmetries for (2^ such that 

pr Yp{Ij) = for all p = 1, . . . ,r and j = 1, . . . , k, 

then Yi, . . . ,Yr are linearly independent generalized symmetries for i28\) . and 

prYp{Ij)=0 for all p=l,...,r and j = l,...,k; 

iv) if under the assumptions of (ii) and (Hi) the symmetries Yi, . . . ,Ys, where s <r, span an involutive distribu- 

s 

tion, i.e., [Yp, Yq] = ^ Cpq{ai, . . . , Ofc, /i, . . . , /fc, Ji, . . . , Jm)yg for all p, q = 1, . . . , s (the commutator is now given by 

9=1 

(3D) !), then the symmetries Yi, . . . ,Ys also span an involutive distribution, [Yp, Yg] = YI '%q^g; foi" all p,q = 1, . . . , s, 

9=1 

where c^q = Cpg(/i, . . . , Ik,bi, . . . ,bk, Ji, ■ ■ ■ , Jm)- If ^ ^ o,'^^ constants (in particular, they do not depend on aj, 
j = 1, . . . ,k), and thus Yg, g = 1, . . . , s, form a Lie algebra, then Yg, g = 1, . . . , s, form an isomorphic Lie algebra. 

Example 4. Consider equation 6.45 from [9j 

d^u fduV 

-d^-=\Tz) 

which admits an integral of the form 

/ = {{du/dzf + (a(l + 2cu))/{2c^)) exp(-2cn). 

Let ai = a and bi = b. Then we have 

/ = 2c2(6exp(-2cn) - {du/dzf + 2cu), 

which is an integral for the transformed equation, 

d'^u/dz^ = c{du/dzf{l - 2c^/{l + 2cu)) + 260^ exp(-2cti)/(l + 2cu). 

Just as for dynamical systems ([T]), for ODEs we also often can add parameters by hand through changes of 
variables, and apply the generalized Stackel transform to the resulting equations. 
Example 5. Consider equation 7.7 from [9], 

d^u 1 d'^u du 2 
dz^ u dz^ dz dz ' 
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which admits an integral of the form 

and rescale u ^ au. This yields the equation 

d^u 1 dP'u du 2 du 
dz^ u dz^ dz dz 

with an integral 

1 d'^u a?v^ 

Again let ai = a and hi = h. Then we have 

which is an integral for the transformed equation, 

dz^ \ u dz"^ \ \ u dz"^ J J j dz 

7 Conclusions and discussion 

In this paper we extend the multiparameter generalized Stackel transform, or the coupling constant metamorphosis, 
to general dynamical systems ([1]) and ODEs and studied the properties of this extension. In particular, we present 
sufficient conditions under which the transformed system inherits the integrability properties of the original one: the 
existence of Lax representation, (sufficiently many) integrals of motion and symmetries, etc. In contrast with the 
Hamiltonian case [T7j, for general dynamical systems ([T]) we can avoid introducing the reciprocal transformation for 
(the solutions of) the equations of motion. As a result, the relationship among the solutions of the original system 
([1]) and the transformed system dH) is much simpler than in the Hamiltonian case. Note that the same approach was 
successfully applied to discrete dynamical systems, see [16] for details. 

Our results naturally lead to a number of open problems related to the generalized Stackel transform, of which 
we list below just a few. 

First of all, it would be very interesting to find out (both in the Hamiltonian and the non-Hamiltonian case) when 
the transformed dynamical system is algebraically [191 [20] completely integrable provided so is the original system. 
On a related note, the study of relationship among the differential Galois groups of the variational equations (see 
e.g. |18[ El |T2] and references therein for the relevant definitions) for original and transformed systems would be of 
interest too. 

Second, we have just barely scratched the surface by noticing in Proposition [1] that the transformed system 
inherits the existence of a Lax representation from the original system, and e.g. understanding what is the precise 
relationship among the Darboux [T3] and Backlund |10^ ITT] transformations for the original and transformed system 
would certainly be worth the while. 

Third, it is highly desirable to study more systematically the issue of when inserting the parameters 'by hand' 
(cf. the above Examples 2 and 5) and subsequent transforming of the resulting systems leads to interesting new 
examples. We expect this technique to yield a significant extension of the pool of exactly solvable dynamical systems 
and ODEs. A good starting point here could be e.g. to find the transformed counterparts for the ODEs that linearize 
on differentiation [6]. 

Finally, the simplest but perhaps also the most important (cf. e.g. [T7] for the Hamiltonian systems) special case 
when integrals of motion are linear in the parameters undoubtedly deserves to be explored in far more details. 
We hope that the present paper will stimulate further research in these and related areas. 



1 d'^u 

ud^^Y' 
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